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Abstract Cross-laminated timber (CLT) panels consist of 
several layers of lumber stacked crosswise and glued 
together on their faces. Prototype sugi CLT floor panels 
were manufactured and bending tests were carried out 
under the different parameters of lumber modulus of 
elasticity (MOE), number of layers, thickness of lumber 
and thickness of CLT panels. On the basis of above tests, 
bending stiffness and moment carrying capacity were 
predicted by Monte Carlo method. MOE of lumber was 
measured by using grading machine and tensile strength of 
lumber was assumed to be 60 % of bending strength based 
on the obtained bending test. Bending stiffness El of CLT 
panels could be estimated by adopting composite theory 
and equivalent section area. Experimental moment carry¬ 
ing capacity showed 12 % higher value than the calculated 
moment carrying capacity by average lumber failure 
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method, and also showed 45 % higher value than the cal¬ 
culated moment carrying capacity by minimum lumber 
failure method due to the reinforcement of the outer layer 
by the neighboring cross layer. 

Keywords Cross-laminated timber • Composite 
theory • Equivalent section area • Moment carrying 
capacity • Monte Carlo method 

Introduction 

Cross-laminated timber (CLT) is produced from lumber 
strips that are stacked crosswise on top each other and 
glued with pressing to form large solid timber elements. 
CLT was developed in Europe few decades ago and has 
been growing up as industrialized engineering wood 
materials. Recently, CLT timber structure has been applied 
to large-scale multi-story buildings in Europe, which 
requires higher seismic and fire-resistant performance. 
Having not yet developed in Japan, CLT structure system 
might be one of the most effective solutions for the use of 
sugi (Cryptomeria japonica ) from the view point of C0 2 
reduction and sustainable forest management. As CLT 
panel is generally used for the floor panel without using 
horizontal members such as beams and joists, the predic¬ 
tion of bending stiffness and moment carrying capacity is 
of great importance in the design of floor elements. 

In this research, prototype sugi CLT floor panels were 
manufactured and bending tests were carried out under the 
different parameters of lumber modulus of elasticity 
(MOE), number of layers, thickness of lumber and thick¬ 
ness of CLT panels. On the basis of above tests, bending 
stiffness and moment carrying capacity were predicted by 
Monte Carlo method. 
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Fig. 1 Typical CLT cross-section of lumber arrangement (W width of CLT panels, T thickness of CLT panel, n and mx number of lumber in 
Y and X directions) 


Predicting bending stiffness and moment carrying 
capacity of CLT panels 

Modeling of lumber arrangement of CLT cross-section 

Figure 1 shows the typical CLT cross-section of lumber 
arrangement. The variables W, T and L are defined CLT 
width, thickness and length as X-, Y- and Z-axis, respec¬ 
tively. The variable n is the number of layers in y-direction 
and variable mx and mz are the number of lumber in X- and 
Z-direction. The variable t t is the lumber thickness of ith 
layer and the E t is the average MOE of lumber of ith layer. 
Main grain direction of lumber at the outer layer is parallel 
to the Z axes. Only the layers having the grain direction 
parallel to that of outer layers are assumed to be effective. 


Equivalent width of each layer and neutral axis 

Equivalent width W cq . k of kth layer is calculated in Eq. (1). 

E k 


Weq-* = W k X 


Eref. 


(i) 


where k is the number of effective layers [1,2] whose grain 
is parallel to the longitudinal direction; W eq . k , the equiva¬ 
lent width of the kth layer; W k , gross width of the kth layer; 
E k , average MOE of the kth layer; and E ref , average MOE 
of reference layer. 


Distance from X-axis to neutral axis is calculated in 
Eq. (2). 

ES (Weq-* x 4) x y, 


N = 


E5 (Weq-* X 4) 


( 2 ) 


where N is the distance from X-axis to neutral axis; Y k , 
distance from X-axis to the centroid of the kth layer; and 
n e ff, number of the effective layer whose grain is parallel to 
the longitudinal direction. 

Calculation of bending stiffness of CLT panel 
by composite theory 

4-nn shows moment of inertia of kth layer calculated by 
the composite theory (parallel layer theory) where the 
longitudinal direction of lumber [3, 4] in Eq. (3). Bending 
stiffness of CLT panel is calculated in Eq. (4). 


4-nn — 4 + A k x N k 


(3) 


where 4-nn is the moment of inertia of the kth layer area 
from the neutral axes; I k , distance from neutral axes to the 
centroid of each kth layer (N k = N-Y k ); and A k , equivalent 
section area multiplied W eq _ k and t k of the kth layer 


«eff 


EIclt — 


Ere f. 


X I t _ 


k- NN: 


( 4 ) 


k= 1 


where EI CLT is the bending stiffness of CLT panel. 
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Calculation of moment carrying capacity 


Moment carrying capacity of glulam is estimated by using 
MOE and MOR properties of lamina that compose the 
glulam [5-8]. Assuming that the layers are not glued each 
other and there is no transmission of shear forces between 
the layers, the bending curvature p of CLT panel shall be 
equal to that of each layer. Radius of curve Up and the 
moment of each layer are expressed in Eqs. (5) and (6). 

1 = _M_ = Mk_ (5) 

p El clj Ekh 


where p is the bending curvature (1/p: radius of curve); M, 
moment of CLT panel; and M k , moment of the Mi layer. 



Evh 


k L k 


El 


M 


CLT 



when the shear stress between each effective layer is 
transmitted by glue, bending stress of each layer is 
expressed in Eq. (7). 


Mk h E k 

CFh k = — = — x - 

Z* Z k EIclt 


l\/T ^k Ek 
x M = — x- 

2 EIclt 


x M. 


(7) 


where a hk is the bending stress of the Mi layer, and Z k , 
section modulus of the Mi layer. 

Tensile stress of each effective layer is expressed in 
Eq. (8). 


&tk — &k * E k 



N k E k 

EIclt 



where o ik , is the tensile stress of the Mi layer and s k , strain 
due to axial force of the Mi layer 

Criterion for failure of CLT panel determined by the 
combined stress of bending and tension of lumber on each 
layer is expressed in Eq. (9). Moment carrying capacity 
M max is calculated in Eq. (10) which is obtained by 
substituting Eqs. (7) and (8) into (9). 


&bk &tk _ j 


fbk ftk 


^nriax — 


El 


CLT 


X 


2 fbk x ftk 


E k (t k x f tk + 2 N k x f hk ) ’ 


(9) 


( 10 ) 


where f bk is the bending strength of lumber at the Mi layer; 
and f tk , tensile strength of lumber at the Mi layer 


Simulation of bending stiffness and moment carrying 
capacity of CLT panels 



Fig. 2 Distribution of lumber MOE by grading machine 



MOE(GPa) 


Fig. 3 Relationship between MOE and MOR of lumber by grading 
machine ( closed circle with finger joint, open circle without finger 
joint). Dashed line 18 % COV with finger joint, dot and dashed line 
18 % COV without finger joint, continuous line average value of 
lamina on JAS glulam, dotted line lower value of lamina on JAS 
glulam 


MOE, bending and tensile strength of lumber 

Figure 2 shows MOE distribution of sugi lumber used for 
CLT panels obtained by continuous grading machine (IIDA 


Kogyo Co., Ltd., Type MGFE-251). Average value of 
MOE along the longitudinal direction of lumber was 
obtained. MOE distribution was fitted to the two distribu¬ 
tion curves, one is the normal probability density 
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Fig. 4 Lumber arrangement of CLT cross-section 
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distribution curve and the other is log-normal probability 
density distribution curve based on the average value and 
standard deviation of 5973 strips of lumber. Normal dis¬ 
tribution fitted to the distribution of lumber better than log¬ 
normal distribution. After grading, lumber was categorized 
to three groups according to the machine grading. Groups 
Gl, G2 and G3 include MOE between 3.5 and 4.5 kN/ 
mm 2 , 4.5 and 6.5 kN/mm 2 and 6.5 and 8.0 kN/mm 2 , 
respectively. 

Figure 3 shows the relationship between MOE and 
modulus of rapture (MOR) of machine-graded lumber. 
MOE and MOR were determined by the bending test of 
lumber according to the bending method C of structural 
Glulam by Japanese Agricultural Standards [9]. Bending 
test specimen of sugi lumber was carried out on the same 
lot as prototype CLT panel production. 

Relation between MOE and MOR showed positive 
correlation and MOR of lumber with finger joints showed 
lower value than that without finger joints. Lower and 
upper range assuming 18 % of COV was shown with 
dotted line. 


Bending stiffness and moment carrying capacity 
by Monte Carlo method 

Figure 4 shows the configuration of CLT cross-section. 
The configuration of simulated CLT panels matched with 
the manufactured CLT panels. CLT panels consisted of 
1000 mm width, 3000 mm length and 3 types of thickness 
of 90 mm, 120 mm and 150 mm. 90-mm-thickness CLT 
panel had three types of layer arrangement consisting of 
three, four or five layers. 120-mm-thickness CLT panel had 
two types of layer arrangement consisting of four or five 
layers. 150-mm-thickness CLT panel consisted of five 
layers. 

The layer including lumber parallel to the longitudinal 
direction had ten pieces of lumber in X-direction and the 
layer including of lumber perpendicular to the longitudinal 
direction had 28 pieces of lumber in Z-direction. 

In E7 type of CLT panel, G3 grade of lumber was used 
for outer layer and Gl and G2 grades of lumber for other 
layers. In Rd type of CLT panel, random Gl, G2 and G3 
grades of lumber were used for each of the layers. 
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Rd type E7 type 




3 4 5 6 7 8 



MOE(GPa) 


Fig. 5 Results of simulation on bending stiffness of CLT panel {open 
circle E7 type, open square Rd type, dotted line MOE of lumber) 

Bending stiffness of CLT panels was simulated by 
Eq. (4) considering variation of lumber MOE. And moment 
carrying capacity was simulated by Eq. (10). 

Determination method of bending stiffness and moment 
carrying capacity of CLT panel 

Bending stiffness and moment carrying capacity of CLT 
panel were calculated as in the following procedure. 

1. Diagram of cumulative distribution curve of normal 
probability density was created by the distribution of 
MOE lumber. 

2. Average and standard deviation values of bending 
strength of lumber was calculated by using the 
relationship between MOE and MOR of lumber with 
finger joints, and diagram of cumulative distribution of 
bending strength of lumber was created from the 
former. 

3. Assuming that tensile strength of lumber is defined by 
60 % of bending strength, a diagram of cumulative 
distribution curve of tensile strength of lumber was 
created. 


With Finger Joint 




A/max (kN* m) 

Fig. 6 Distribution of moment carrying capacity of CLT panel 
calculated by using minimum bending strength of lumber at the outer 
layer of tensile side (open circle El type, open triangle Rd type, 
dotted line 5th percentile lower) 


4. MOE, bending and tensile strength of each lumber was 
chosen by Monte Carlo method using each cumulative 
curve, where configuration of CLT cross-section 
including chosen MOE, bending and tensile strength 
of lumber follows in Fig. 4. 

5. Bending stiffness of the CLT panel was calculated by 
Eq. (4) of which layer’s MOE was the average MOE of 
lumber in that layer. 

6. Moment carrying capacity of CLT panel was calcu¬ 
lated by Eq. (10) of which the outer lumber strength at 
the tension side was decided by the minimum bending 
strength among that of the particular layer (minimum 
strength model) or the average bending strength of the 
concerned layer (average strength model). 

Results of simulation on bending stiffness and moment 
carrying capacity 

Figure 5 shows the results of simulation on bending stiff¬ 
ness of CLT with 90, 120 and 150 mm thickness. Variation 
of MOE of each CLT panel showed smaller value than the 
variation of MOE of lumber consisted of CLT. E7 type 
CLT panel showed higher MOE than the Rd type CLT 
panel. Bending stiffness of CLT panel was affected by the 
MOE of lumber at the outer layer. 
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Figure 6 shows the distribution of moment carrying 
capacity of CLT panels with 90, 120 and 150 mm thick¬ 
ness. Moment carrying capacity was calculated by the 
minimum bending and tensile strengths among the lumber 
in the outer layer (minimum strength model). 

Upper graph shows the results of simulation by the outer 
layer with finger joints and lower graph shows the results 
of simulation by the outer layer without finger joints. 
Maximum moment carrying capacity of CLT panel with 
finger joints showed lower value than the average moment 
carrying capacity of CLT without finger joints. However, 
distribution range of moment carrying capacity with finger 
joints was narrower than that without finger joints. 

Bending test 

Test specimen 

Specimen of bending test was manufactured in the similar 
procedure as the specification of simulation for 1000 mm 
width, 3000 mm length and 90, 120 and 150 mm thickness. 
MOE of the lumber was measured by machine grading as 
shown in Fig. 2. Lumber arrangement of CLT cross-section 
was shown in Fig. 4. Water-based polymer isocyanate 
adhesive was used for lumber lamination and finger join¬ 
ting (15 mm in length, 3.8 mm in pitch, 1/12 in slope, 
0.7 mm in tip thickness. Orientation of finger joints of 
lumber was vertical). And no glue was applied to the edge 
joint of lumber on CLT manufactured. Density and mois¬ 
ture content of CLT panels are shown in Table 1. Moisture 
content of CLT panel was measured by wood moisture 
tester which uses the dielectric constant measuring method 
(Kett Electric Laboratory HM-530). Two specimens were 
tested for 150 and 120 mm thickness of CLT panel and 
only one specimen was tested for 90 mm thickness and the 
90-mm-thickness CLT test specimen was measured only 
for stiffness of panels without failure. 

Test method 

Bending test was carried out according to Japan Indus¬ 
trial Standard JIS A 1414-2:2010 [10]. Bending test 
facility had a 500-kN loading capacity (Maekawa Testing 
Machine Mfg. Co., Ltd. Type IPU-20/100B-B1) located 
in Hiroshima Prefectural Technology Research Institute 
Forestry Research Center. Loads were applied at one- 
third point of 2700 mm span as shown in Fig. 7. Load 
was measured by load cell at the center of loading beam 
and vertical deflections were measured at the center and 
the supporting points. (Tokyo Sokki Kenkyujo Co., Ltd, 
SDP-100C). The average displacement at the supports 
was subtracted from the central displacement. The 


embedding displacements of CLT panels at the supports 
were nevertheless small. Deflection at the moment con¬ 
stant area was measured at the both sides of yoke 
equipment (Tokyo Sokki Kenkyujo Co., Ltd, CDP-50). 
Strains of outer layer lumber on tension and compression 
side were measured by strain gauge (Tokyo Sokki 
Kenkyujo Co., Ltd PFL-20-11-5L). Measuring position 
of strain gauge was shown in Fig. 8. 

Loading protocol at the center deflection of span con¬ 
sisted of one-way cyclic loading of 1/300, 1/150, 1/100, 
1/75, 1/60 and 1/50 of total span 2700 mm. Loading of 
90-mm-thickness CLT panel was stopped at the deflection 
of 1/100 without failure. 

Bending stiffness based on the deflection, bending 
stiffness based on moment constant deflection by yoke 
equipment and maximum moment were calculated by Eqs. 
(11), (12) and (13), respectively. 


23Z/ A P 

El — x 

1296 A (Renter 

(ii) 

L] A P 

Elyoke — X 

16 Adyoke 

(12) 

P max X L 

M — 

1V1 max — s 

6 

(13) 


Table 1 Density and moisture content of CLT panels 


Thickness 

(mm) 

Layer 

Arrangement 
of lumber 

Number 

of 

specimen 

Moisture 

content 

(%) 

Density 

(kg/m 3 ) 

150 

5 

Rd 

2 

9.1 

436 





12.2 

431 



E7 

2 

12.4 

442 





15.0 

431 

120 

4 

Rd 

2 

11.3 

442 





13.3 

439 



E7 

2 

10.7 

433 





12.8 

436 


5 

Rd 

2 

11.4 

444 





10.8 

439 



E7 

2 

12.4 

436 





10.8 

442 

90 

3 

Rd 

1 

12.7 

441 



E7 

1 

13.0 

433 


4 

Rd 

1 

15.3 

441 



E7 

1 

8.1 

448 


5 

Rd 

1 

12.6 

448 



E7 

1 

12.8 

441 




Average 

12.0 

439 




Standard 

deviation 

1.8 

5 
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Fig. 7 Setup of the test 
specimen for bending (unit, 
mm) (setup the hard wood 
plate 160 x 1300 mm at 
loading point and steel 
plate 200 x 1300 mm at 
supporting point.) 



Fig. 8 Measuring position of 
strain by strain gauge (unit, 
mm) 
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where El, bending stiffness on deflection by 2.7 m span 
(kN m ); EI yoke , bending stiffness on yoke deflection 
(kN m ); A P, AS center and A^ yoke are load and elastic 
range of deflection difference at center and yoke; L, span 
of bending (L = 2.7 m); L 1? span of bending between the 
moment constant area (L\ = 0.9 m); M max , moment 
carrying capacity (kN m); and P max , maximum Load 
(kN). 


Results of bending tests 

Figure 9 shows the relationship between the moment and 
deflection of CLT panel consisted of 120 mm (5ply), 
120 mm (4ply) and 150 mm (5ply). Stiffness and strength 
of E7 type showed higher value than the Rd type. Test 
results of two specimens showed almost the similar per¬ 
formance. Under cyclic loading below L/60 (center 
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T120-5Ply 


T120-4Ply 


T150-5Ply 





Deflection (mm) Deflection (mm) 
Fig. 9 Relationship between moment and deflection of bending test 


CLT 120mm 5Ply Type-E7 



Position of Strain Gage 


CLT 150mm 5Ply Type-Rd 
Compression side 


X X x 



§ s 




SGI SG2 SG3 SG4 SG5 
Position of Strain Gage 


0:1/300, 0:1/150, Ckl/lOO, A:l/75, V:l/60, (Q): 1/50 

Fig. 10 Typical strain distribution of CLT panel under bending test. 
Continuous line E7 type, dotted line Rd type 


deflection of 45 mm), reduction of stiffness of CLT panels 
was not observed. Typical strain distribution of CLT panel 
under bending test was shown in Fig. 10. Left figure shows 
E7 type CLT panel (120 mm and 5ply) and right figure 
shows Rd type CLT panel (150 mm and 5ply). Within the 
small deflection of L/60, strain distribution in X-direction 
of outer layer was almost constant on both tension and 
compression sides of CLT panel. When finger joint of 
certain lumber failed on the tension side of outer layer, the 
strains of adjacent lumber did not decrease. Another lum¬ 
ber of tension side next to the failed lumber shared the 
strain. 

Figure 11 shows relationship between bending stiffness 
calculated by center deflection of 2700 mm span and cal¬ 
culated by yoke deflection. Both bending stiffness shows 
high correlation. Bending stiffness calculated by yoke 
deflection showed 8 % higher than the bending stiffness 



3 2 

Exprimental El on Center Deflection(x 10 kNm~) 

Fig. 11 Relationship between the bending stiffness based on the 
deflection of 2700 mm span and El based on the yoke deflection 

calculated by center deflection on 2700 mm span. It is 
because bending stiffness calculated by center deflection of 
2700 mm span included not only bending, but also shear 
deflection. 

Figure 12 shows typical failure of CLT panel on the 
tension side. At first, one of the finger joint at the tension 
side of constant moment area failed and then the next 
lumber failed due to the defects such as the finger joints, 
knot of lumber and cross grain continuously. 

Comparison of simulation and test results 

Figure 13 shows the relationship between the calculated 
and experimental bending stiffness. Number of calculation 
by Monte Carlo method was 300 on each CLT panels. 
Therefore, each plot shows average of calculated bending 
stiffness and horizontal line shows 3 a ( a , standard devia¬ 
tion). Circle plot and square plot shows experimental 
bending stiffness calculated by yoke deflection and by 
center deflection on 2700 mm span. Calculated bending 
stiffness based on the lumber MOE by grading machine 
with composite theory and equivalent section area was 
concordant with experimental bending stiffness. Experi¬ 
mental bending stiffness by measuring yoke deflection 
fitted more than the experimental bending stiffness by 
measuring center deflection by 2700 mm span. 

Figure 14 shows the relationship between the calculated 
and experimental moment carrying capacity. Circle plot 
shows calculated moment used average strength of lumber 
at the outer layer on the tension side and square plot shows 
calculated moment used minimum strength of lumber at the 
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Fig. 12 Typical failure of CLT panel on tension side (type Rd 
150 mm 5ply) 



Fig. 13 Relationship between calculated and experimental El on 
CLT panel 

outer layer on the tension side. Experimental moment 
carrying capacity showed 12 % higher value than the cal¬ 
culated moment carrying capacity by average lumber 
failure method and showed 45 % higher value than the 
calculated moment carrying capacity by minimum lumber 
failure method. Those results indicate that neighboring 
cross layer distribute bending and tensile stress over the 
outer layer. 

Precise prediction of moment carrying capacity of CLT 
panel by Monte Carlo method is effective. However, in the 
simulation, it is necessary to prepare the accurate lumber 
strength based on the MOE by grading machine. Assuming 
that the CLT is composed of all the minimum MOE by 



Fig. 14 Relationship between calculated and experimental maximum 
moment on CLT panel 



Mmax by DDM (kN* m) 

Fig. 15 Relationship between the moment carrying capacity by 
DDM and moment carrying capacity of 5th percentile lower by Monte 
Carlo method 

grading, it is possible to simplify to design CLT panels. 
Deterministic design method (DDM) deals with the model 
whose MOE and MOR of all lumber in each layer is iden¬ 
tical and based on the 5th percentile lower limit of lumber. 
MOE of 6.0 GPa and 22.5 MPa and 13.5 MPa for bending 
and tensile strength were used for E7 type of CLT panel. 
MOE of 3.5 GPa and 17.0 MPa and 10.2 MPa for bending 
and tensile strength were used for RD type of CLT panel. 
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Figure 15 shows the relationship between the moment 
carrying capacity by deterministic design method (DDM) 
assuming that the CLT is composed of all the minimum 
MOE lumber and moment carrying capacity of 5th per¬ 
centile lower by Monte Carlo method. Fifth percentile 
moment carrying capacity by minimum strength model 
coincided well with the moment carrying capacity calcu¬ 
lated by DDM. Fifth percentile moment carrying capacity 
by average strength model showed 1.6 times higher than 
both DDM and minimum strength model. Considering the 
simulated moment carrying capacity by average strength 
model was close to the experimental values in Fig. 14, 
DDM and minimum strength model underestimated the 
moment carrying capacity of CLT panel. 

Conclusions 

The following results were obtained in this study. 

1. Bending stiffness calculated with composite theory 
and equivalent section area by Monte Carlo method 
was concordant with experimental bending stiffness. 

2. Experimental moment carrying capacity showed 12 % 
higher value than the moment carrying capacity 
calculated by average strength model and showed 
45 % higher value than the moment carrying capacity 
calculated by minimum strength model. It is supposed 
that the cross layer connecting each lumber of outer 
layer transmits the bending and tensile stress and 
contributes to increase the moment carrying capacity 
of CLT panels. 

3. Within the small deflection of L/60, strain distribution 
in the lateral direction of outer layer was almost 
constant on both tension and compression sides of 
CLT panel. The fact that the strains of adjacent lumber 
to the failed lumber with finger joint did not decrease 
shows that the adjacent lumber shares the tensile and 
bending stress because of the existence of the cross 
layer. 

4. Moment carrying capacity of CLT panel obtained from 
the 5th percentile lower limit calculated with minimum 
strength model was concordant with that calculated 
with DDM, and that obtained from the 5th percentile 
lower limit calculated with average strength model 
showed 1.6 times higher value than that calculated 


with DDM which was close to the experimental 
results. This indicates that the moment carrying 
capacity calculated by DDM underestimates the 
moment carrying capacity of CLT panel although it 
gives safety side design values. 
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